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■ Abstract 

o : 

£SJ , We consider f(R,T) theory of gravity, where R is the curvature scalar and T the trace of the 

energy momentum tensor. Attention is attached to the special case, f(R,T) = R + 2/(T) and two 
expressions are assumed for the function f(T), ^„|^ and 03 In 9 (&3T" 1 ), where 01, 0,2, bi, 62, n, 
. 03, 63, q and m are input parameters. We observe that by adjusting suitably these input parameters, 

t-H ■ energy conditions can be satisfied. Moreover, an analyse of the perturbations and stabilities of de 

1— 1 ■ Sitter solutions and power-law solutions is performed with the use of the two models. The results 

^jH- show that for some values of the input parameters, for which energy conditions are satisfied, de Sitter 

' solutions and power-law solutions may be stables. 

Pacs numbers: 04.50.-h, 04.50.Kd, 98.80.-k 

^t- : 

> . 
00 ■ 

t-» ■ 1 Introduction 

: 

It is well known that General Relativity (GR) based on the Einstein-Hilbert action (without taking 
into account the dark energy) can not explain the acceleration of the early and late universe. Therefore, 

<N ■ 

^-H ' GR does not describe precisely gravity and it is quite reasonable to modify it in order to get theories that 

admit inflation and imitate the dark energy. The first tentative in this way is substituting Einstein-Hilbert 
term by an arbitrary function of the curvature scalar R, this is the so-called i(R) theory of gravity. This 

5_i , 

theory has been widely studied and interesting results have been found PQ-[2]. In the same way, other 
alternative theory of modified gravity has been introduced, the so-called Gauss-Bonnet gravity, /(G), as 
a general function of the Gauss-Bonnet invariant term G [3j . Other combinations of scalars are also used 
as the generalised f(R,G) and f(R,P,Q) [IE], where P = R^W and Q = R t2lJrTT R^ CTT (here R^ 
and R^mt arc the Ricci tensor the Riemann tensor, respectively). 

In this present paper, attention is attached to a type of the so-called f(R, T) theory of gravity, where 
T denotes the trace of the energy momentum tensor. This generalization of f{R) gravity has been made 
first by Harko et al [6j . The dependence on T can take source from the introduction of exotic imperfect 
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fluid or from quantum effect (conformal anomaly). In this theory, the equations of motion show the 
presence of an extra-force acting on the test particles, and the motion are generally non-geodesic. This 
theory also relates the cosmic acceleration, not only due to the contribution of geometrical terms, but 
also to the matter contents [6] . Some results have been obtained with this theory. In [7] , the cosmological 
reconstruction of f(R,T) describing transition from matter dominated phase to the late accelerated 
epoch of the universe is performed. Also in the same way for exploring cosmological scenarios based on 
this theory, f(R,T) function has been numerically reproduced according to holographic dark energy |S]. 
Moreover it is shown that dust reproduces ACDM , phantom-non-phantom and the phantom cosmology 
with f(R,T) theory j9]. The general technique for performing this reproduction of AC DAI model in 
FKW[f| metric cosmological evolution is widely developed in [U [TU]. The f{R,T) models that are able 
to reproduce the fourth known types of future finite-time singularities have been investigated [11] . The 
authors also introduced conformal anomaly as quantum effects near the singularities and observe that it 
cannot remove the sudden singularity or the type IV one, but, for some values of a input parameter, the 
big rip and the type III singularity may be avoided from the effects of the conformal anomaly. Moreover, 
they found as interesting result that, even without taking into account conformal anomaly, singularities 
(the Big rip and the type-Ill) may still be removed thanks to the presence of the T contribution of the 
f(R,T) theory. 

Note that singularities appear when energy conditions are violated. Our task in this paper is to check 
the viability of some models of f(R,T) according to the energy conditions. The energy conditions arc 
formulated by the use of the Raychaudhuri equation for expansion and is based on the attractive character 
of the gravity. We refer the readers to Refs [12]- [17]. where energy conditions are widely analysed for the 
cosmology settings, in f(R) and f(G) gravities. 

In this paper, we assume a special form of f(R, T), that is, f(R, T) — R + 2/(T), the usual Einstein- 
Hilbert term plus a T dependent function f(T). Two expressions of f(T), and a 3 In 9 (b 3 T m ) 

are investigated, where a%, a,2, b%, 62, n, a 3 , 63, q and m are parameters to be suitably adjusted in order 
to obtain viable models according to energy conditions. The first expression is similar to that used in 
[E1QJ5], where, at the place of the trace T, the Gauss-Bonnet invariant is used. The main motivation of 
this choice is explained in the section 4. The second form is chosen due to its interesting aspect, because 
its corresponding f(R,T) model avoids itself the appearance of the type-I finite-time future singularity 
(Big Rip) [11] . thanks to the ordinary trace T-terms contributions. Thus it is quite reasonable to check 
for which values of the input parameters this model can be acceptable as cosmological model. 

In order to reach the acceptable cosmological models, we analyse the perturbations and stabilities of 
de Sitter solutions and power-laws solutions in the framework of the special R + 2/(T) gravity, by using 
the two models proposed in this work. We observe that for some values of the input parameters, for both 
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models, the stabilities of de Sitter solutions and power-law solutions are realized and compatibles with 
some energy conditions and the late time acceleration of the universe. 

The paper is outlined as follows. In section 2, we briefly present the general formalism of the theory, 
putting out the general equations of motion for a f(R,T) = fi(R) + f2(T) gravity, where fi(R) and 
f2{T) are respectively function of the curvature scalar and the trace of the energy momentum tensor. 
The section 3 is devoted to the general aspects of the energy conditions. The f(R,T) = R + 2/(T) 
gravity is assumed in the section 4, where the two functions considered for f(T) are studied, putting out 
the conditions on the input parameters for obtaining some viable models of f(R,T). The perturbations 
and stabilities of de Sitter and power-law solutions are investigated in the sections 5. Discussions and 
perspectives are presented in the section 6. 



2 General formalism 

Let us assume the modified gravity replacing the Ricci scalar R in Einstein gravity by an arbitrary 
function /(i?,T), and writing the total action as 

S= ^2 J dix V^g[f(R,T) + C] , (l) 

where k 2 = 8irG, G being the gravitational constant and T = g^ v T^ v the trace of the matter energy 
momentum tensor which is defined by 

2 SU=gC) 

T — — - - (21 

This modified gravity theory has been considered first in [B] and the equations of motion, using the metric 
formalism, have been explicitly obtained as 

f R (R, T)R MV - \f{R, T)g MV + (g^D - V M V W ) f R (R, T) = k 2 T^ - f T (R, T)T^ - f T (R, T)Q MV , (3) 

where fu(R,T) and fr{R,T) denote the derivative of f(R, T) with respect to R and T, respectively, and 
O^u is defined by 

A " y ~- 9 Sg^ ~ ^J-^ + 9^L, zg dg ^ dgaf} ■ W 
Let us consider the stress-energy tensor of the matter as 

Tfj, v = (p + Pju^Uv - pg^ . (5) 

This expression is a consequence of a statement of the fluid mechanics for which a possible way for 
treating the scalar invariant C is to interpret it as the pressure p of the fluid, i.e. C = —p (the negative 
sign is due to the signature used here) UM ■ Here, p, p and u M are respectively the energy density, the 
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pressure and the four- velocity of a perfect fluid considered as the classical matter content of the universe. 
Hence, Eq.(j4|) becomes 

6^ = -2T MI/ - pg^ . (6) 

From now, we will treat the models of type f(R, T) = fi(R) + f2(T) and set k 2 = 1. Then, the equation 
([3]) becomes 



fi R (R)R»v - \h{R)g„v + (g^O - f 1R (R) = + f 2T {T)T^ 



/»(!>+ ~/ 2 (T) 



6V ( 7 ) 



3 Energy conditions 

The energy conditions are essentially based on the Raychaudhuri equation that describes the behaviour 
of a congruence of timelike, spacelike or lightlike curves. It is commonly used to study and establish 
singularities of the spacetime. Considering a congruence of the curves fulfilling the spacetime manifold 
and also that there are expansion, distortion and relative rotation between the curves of the congruence, 
then the spacetime manifold described by the congruence is curved. For the purposes of this work we will 
just consider the timelike and spacelike curves for which the Raychaudhuri equation reads, respectively 

[201 EO 

R» V V»V V + + - uvu^ + = , (8) 



1 d6 

R^k v + -9 2 + o^ v - + — = , (9) 



where 6 is the expansion scalar describing the expansion of volume, r and A are positive parameters 
used to describe the curved of the congruence, a^ v the shear tensor which measures the distortion of 
the volume, uj^ v the vorticity tensor which measures the rotation of the curves, and and are 
respectively timelike and lightlike vectors tangent to the curves. In this work, we are interested to the 
situation for small distortions of the volume, without rotation, in such a way that the quadratic terms 
in the Raychaudhuri equation may be disregarded (they are like second order corrections). Then, the 
equation can be integrated given the scalar of expansion as a function of the Ricci tensor: 

9 = -t^FF = -XR^k" . (10) 

The condition for attractive gravity is 9 < 0, imposing R 11U V^V 1 ' > and R ilv k^k u > 0. These two 
conditions are called the strong and null energy conditions, respectively. 

For equivalence to GR, by just dividing by fi(R) (different from zero), one can cast Eq. ([7]) in the 
following form 

R^-\Rg^ = T;l f , (li) 
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where the effective energy momentum tensor T^J is defined by 

T e JJ = + fMT)T^ + \ [2f 2T (T)p + f 2 (T) + MR) - Rfm(R)} 

-( g ^n-v^v„)f 1R (R)} ■ (12) 

Since the vector field k^ is lightlike, one gets k^k^ — 0. Then, by multiplying Eq. (TTTj) by the bi-vector 
k^k v , the null energy condition comes to be equivalent to T^Jk^k v > 0. Remark also that Eq. (fTTj) is 
equivalent to 

R»u = T*JJ - l -T e ffg^ v . (13) 

Then, the strong energy condition reads [TJjff — \T e * ' g^ V^V V > 0. By assuming that the total 
content of the universe behaves as perfect fluid, Eq. ([5]) also holds, by just substituting p and p, by p e // 
and Peff respectively (the effective energy density and effective pressure). Thus, the null energy condition 
imposes 

TtffWW = {peff + Peff) U fi U u k"k V - Peffg^k" > . (14) 

Since, k^ is a lightlike vector, one has g^k^k 1 ' — k^k^ = 0. Also, we always get u tl u u k^k 1 ' — (u^k 1 - 1 ) 2 > 0. 
Thus, the null energy condition for the effective perfect fluid reduces to 



Peff +Peff > . (15) 



For the strong energy condition, one has 



- -T 9la> ) rr 



(Peff + Peff) U^U V - Peffg^v ~ ^ (Peff ~ iPef f) 9^ 



VV^ 1 > ,(16) 

= ( P eff+5peff)V^ (17) 

where we made use of T e ^ = p e ff — Speff and u tl u iy V IJ 'V v = V^V^ for obtaining the second equality|f| 
Remark that the condition (TTE| is equivalent to 



(Peff + Peff) V"% > \ (Peff ~ Peff) V% , (18) 

which imposes (p e ff + Pef f)V^V^ > 0. Since the vector is timelike, one has V^V^ > 0, and the 
strong energy condition reduces to 

Peff +Peff > , Peff + Zpeff > . (19) 

Note that the condition (Tffi - \T e ^ g^) V»V > is equivalent to T e JJV^V v > ±T eff V»V v . This 
requires to have the necessary condition T^V^V" > 0. The weak energy condition requires to this 



6 Since is an unitary timelike vector (four- velocity) and the metric component goo = 9 00 = 1 ( scc Eq. H24H - ). one gets 
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condition holds for any non-spacelike vector [201 [22]. This means that the condition holds for both 
lightlike and timelike vectors. The use of the lightlike vector is equivalent to the null energy condition as 
previously demonstrated. Now, for a timelike vector, one gets 

T «JJV»V V = [(p eff + Peff ) Ull u v ~ p eff9 ^\ W > , 

= PeffWV,, > . (20) 

Since, is a timelike vector, one always has V^V^ > 0. Thus, the weak energy condition for the effective 
perfect fluid reads 

Peff>0, (Peff+P*ff)>0 ■ (21) 

For the dominant energy condition, the required condition is that, besides the weak energy condition, for 
any timelike vector V^, T^V^ is a non-spacelike vector. In other word, this means that no signal can 
propagate faster than light. It follows that Peff < Peff, since, for any know matter, pressures are small 
when the energy density is small. Thus, the dominant energy condition results in 

Peff > , p e ff + p e ff > , p e ff ~ Peff > . (22) 

Therefore, the energy conditions, as known in GR, can also be applied in this modified theory of 
gravity by substituting the ordinary energy density p and pressure p in GR by the effective ones, p e ff 
and p e ff. 

In what follows, we will consider models of type f(R,T) = R + 2/(T), i.e., the usual Einstein-Hilbert 
term plus trace depending term 2/(T). This amounts to consider fi(R) = R and fziT) — 2/(T). The 
factor 2 is used just for letting the field equations more easier to be treated. We will also assume that 
the ordinary content of the universe is pressureless and satisfies the energy conditions (just p > 0). 



4 Testing some f(R,T) = R + 2/(T) models from energy condi- 
tions 

In this section we will present the conditions required on p and the algebraic function f(T) for realizing 
each type of energy conditions. For this end, we first need to establish the respective expression of the 
effective energy density p e // & n d effective pressure Peff- According to the assumptions made at the end 
of the previous section, Eq. (0 becomes 

Rnu - \R9^ - + 2/ T (T)T M „ + f{T) 9tlv . (23) 

Considering the flat FRW space-time described by the metric 

ds 2 =dt 2 -a 2 (t)dx 2 , (24) 
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where a(t) is the scale factor. The 00 and it components of (|23|) can be written as 

3H 2 = Peff , (25) 

-2H-3H 2 = p eff , (26) 

where the effective energy density and pressure are defined as 

p eff = p + 2pf T (T) + f(T) , (27) 

Peff = -f(T). (28) 

By using the above expressions of the effective energy density and pressure, we get the null energy 
condition (NEC), the weak energy condition (WEC), the strong energy condition (SEC) and the dominant 
energy condition (DEC) by 



NEC: p e ff + peff =p[l + 2f T (T)] >0 , (29) 

WEC: Peff =p + 2pf T (T) + f(T)>0, p e ff+Peff>0, (30) 

SEC: p eff + 3 Peff = p + 2pf T (T)-2f(T)>0, p e// + Pe// >0, (31) 

DEC: p ef f-p eff = p + 2pf T (T) + 2f(T)>0, p e/ / + p e// > , p eff > . (32) 

We propose to test two models of f(T) in the way to make them satisfying the energy conditions. For the 
first model, we consider a function f(T) such that for large and small values of the trace, it converges. 
We start by assuming first this function as f(T) = °^„|^ , where ai, a,2, &i, 62 and n are parameters 
to be adjusted in order to obtain models that satisfy some or all the energy conditions. Note that this 
expression is similar to that used in [18l [19], where, at the place of the trace T, the Gauss-Bonnet 
invariant is used. A motivation in using this model is that when n > 0, it prevents divergence for 
large and small values of the tracj^. Remark that for T — > 00, /(T) = a±/a2 is finite, and for T — > 0, 
f(T) = 61/62 which is also finite. However, when n < 0, the situation inverts, and for large and small 



values of the trace*], one obtains f(T) = 61/62 and f(T) = a\/a2, respectively, and the divergence is 
still prevented. But note that all this is a mathematical point of view, which also may be useful at the 
physical point of view. Let us now explore the cosmological feature of this model. Here we can discuss 
the occurrence of cosmic acceleration, which may impose some constraints to the parameters, reducing 
the degree of freedom of the model. Searching for the power-like solutions, one has: a — aot x , H = it -1 , 



7 This requires to fix ai 7^ for large values of trace and 62 7^ for small values of the trace. 
8 This requires to fix ai ^ for small values of trace and 62 7^ for large values of the trace. 
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H = -xt- 2 and f(T) -> g(t) = (a^ 3 ™ + b 1 )/{a 2 t~ 3xn + b 2 ), with a. t = a z p Q a^ 3 (i = 1,2), where 
we used T = p = pa~ 3 (because the ordinary content of the universe is dust). With the scale factor 
being used here, the early universe corresponds to t — > and the late time correspond to t — > oo, where 
x > 0. At early time, the curvature scalar R — > oo (where the f(T) contribution may be neglected) and 
at late-time, the universe may be characterized by the AC DAI model, i.e., our model must behave like 
R + 2A. For obtaining this feature, for n < 0, the cosmological constant reads A = axja-i — ai/a^, while 
for n > 0, one gets A = &i/&2- By calling w e // the effective parameter of the equation of state, one has 
w e// = Peff/Peff = ~ f (T) / (p + 2pf P + f(T)). The requirement of guaranteeing the acceleration of the 
universe, without falling into phantom model, is —1 < u) e tt < —1/3. Observe that at time, (t — > oo), 
uj e ff = — 1 (this holds only when n > 0). Thus, we see that some values of the input parameters can make 
the model providing the late time acceleration, i.e., for n > 0, ai/tt2 > and &1/&2 = A. It is important 
to note that when the acceleration is guaranteed, at least the strong energy condition is violated. But 
here, we recall that the aim of this work is to find the range of the parameters for which the energy 
conditions are satisfied. We will develop this in the subsection 4.1. 

For the second model, we assume f(T) = 03 In 9 (b^T" 1 ) where 03, m and q are also adjustable param- 
eters and 63 is assumed to be positive and non null. This form is chosen due to its interesting aspect, 
because it is a type of the f(T) model for which, the future finite-time singularity of type-I (the Big Rip) 
can be cured in f(R,T) = R + 2/(T) gravity 11 . Despite the aim of this work not being the study of 
the avoidance of singularities, we can try to show a little more how this type of model can prevent the 
Big Rip. Remember that the Big Rip is the type of the future finite-time singularity for which, the scale 
factor, the energy density and the pressure diverge. This also implies that at the singularity time, here 
denoted by t s , the curvature scalar R diverges. Then, in order to show the avoidance of the Big Rip from 
the contribution of the /(T), let us take the trace of the Eq. (f23|) . i.e. 

-R = p + 2pf T (T)+Af(T) . (33) 

Remark that in |llj . the Hubble parameter is written as H — h(t s — t)~~ 01 , where h is a positive real 
parameter and a > 1; this is the condition for the appearance of the Big Rip. In this case, the scale 
factor reads a = aexp [h(t s -t) 1 - a /(a-l)], where a is a positive parameter. Hence, as the singularity is 
approached, the curvature R behaves like (t s — t)~ 2a . On the other hand, one can determine the quantity 
b 3 T m = b 3 p m = b 3 p m a~ 3m exp [3hm(t s ~ t) 1_a /(l - a)] and the f(T) = f(p) behaves like (t s - 1)^ 1 '^. 
From these points, with m < 0, we see that for q > 2a/ (a — 1), q > and q(l — a) < because of a > 1. 
Moreover, we can see that q(l ~ a) < —2a. Hence, as the singularity is approached, the T contribution 
diverges more than the curvature scalar, and the Big Rip may be avoided, without the need of quantum 
effects. This proves that the model may allow itself the avoidance of the Big Rip. Here, we propose to 
determine the range for the parameters for which the energy conditions are satisfied, and guaranteeing 
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the auto avoidance of the Big Rip. Note also that in [TT] it is demonstrated that this model can provide 
the late time accelerated expansion of the universe. 



4.1 Studying the case f(T) = |^±| 

Our task here is to put out the constraints on the input parameters in order to get a R + 2/(T) type 
model that satisfies the energy conditions. According to the sign of the parameter n, and assuming that 
d2 and 62 cannot be identically null, the model can be cast into two different forms. In fact, for the late 
time stage of the universe, by dividing the parameters of the model by 02 (n > 0) and 62 {n < 0), one 
gets respectively the models f(p) = (Ap n + B 1 )/(p n + B 2 ) and f(p) = {A x p n + A)/(A 2 p n + 1), where 
the cosmological constant is characterized by 01/02 (for n > 0) and 61/62 (for n < 0), and A\ = ai/62, 
A 2 = 02/62, B\ = bi/a 2 and B 2 = b 2 /a 2 . In this case, the model which initially was four parameters 
dependent, under the cosmological constraints, becomes three parameters dependent, A, B\ and B 2 for 
n > 0, and A, A\ and A 2 for n < 0. Since the cosmological constant is known 14^|, the model turns into 
two parameters dependent. 

The first derivative of f(T) with respect to T (or the derivative of f(p) with respect to p) reads 



n{KB 2 -B 1 )p- 



1- 1 



for n > 



• The NEC 

Since we have assumed that the ordinary content of the universe satisfies all the energy conditions, 
the condition (29) reduces to 1 + 2/ T (T) > 0, (or 1/2 + f p (p) > 0). One can calculate 1/2 + f p (p) as 



1 J 2 " P "' 1(A g: + B ff P " +B2)2 , for n > 

I 2(A 2P "+i) a ' toi n < U , 

whose the sign can just be characterised by that of the numerator, since the denominator is always positive. 
If we take the numerator as a function of the ordinary energy density p and the input parameters, we 
just need to analyse the sign of this latter. The evident conditions ^ for which the numerator is positive 
are presented as follows: 

* B x > 0, B 2 > 0, B x /B 2 < A for n > , 

* A x > 0, A 2 > 0, B 1 /B 2 > A for n< , 

* Si < 0, B 2 > 0, for n > , 

* Ai < 0, A 2 > 0, for n < . 



9 The cosmological constant is positive and is A ~ 1.7 X 10 121 Planck units. 

10 Here, we call "evident conditions" the conditions based just on the signs of the parameter without the need of knowing 
their absolute values. 
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Indeed, the above conditions lead to the positivity 2np n ~ 1 (AB 2 —B\) > for n > and 2np n ~ 1 (Ax — 
A 2 A) > for n < 0. Observe that there are still situations in which the above quantities are negative 
but the numerators in (|35|) continuing positive, i.e., 

* A 1 > A 2 A and \2np n - x (Ax - A 2 A) \ < \ (A 2 p n + i f | for n < , 

* AB 2 < B x and \2np n - x (AB 2 — -B x ) | < j (p n + B 2 f | for n > . 

In these cases, one can plot the function in terms of two of the parameters, fixing the other. Despite 
knowing the sign of the considered parameters with what respect the function may be plotted, the 
important here is their rank, i.e. the interval to which they must belong in order to produce the positivity 
of the function. Some examples are presented in Fig. 1. 

* The WEC 

This condition is realized when the NEC is, plus the condition p e ff > 0. Note that the complete 
expression and condition of the NEC read 

2n(AB 2 -B 1 )p' l + p n+1 + 2B 2 p n+1 +B 2 2 p>0 , for n > , (36) 
2n {A x - A 2 A) p n + A 2 2 p n+1 + 2A 2 p n+1 + p > , for n < . (37) 

These expressions are obtained by multiplying the numerators in p5[) by p. We didn't need to use this 
complete expression for determining the conditions on the input parameters in the case of the NEC, since 
the ordinary energy density is assumed as positive quantity. Besides to (|36j) -(|37 |) . the second condition 
for satisfying the WEC is 

Peff = p + 2pf p + f(p) >0 , (38) 

having in mind that the ordinary content is assumed as pressureless. By using f(p) according to the 
functions in (|3"5)) becomes 

Ap 2n + p 2n+1 + 2B 2 p n+1 + (2nAB 2 - 2nB 1 + AB 2 + Bi) p n + B\p + B X B 2 > , for n > 0, (39) 
A x A 2P 2n + A 2 p 2n+1 + 2A 2 p n+1 + (2nA 1 - 2nA 2 A + A x + AA 2 ) p n + p + A > , for n < 0. (40) 

Note here that we just use the numerator of the fractions whose the denominators are always positives. 
By combining (J36J) with ([39]) -g0]), one gets for the WEC 

Ap 2n + 2p 2n+1 + AB 2 p n+1 + 2B 2 p + [(1 + 4n)AB 2 + (1 - 4n)Bi] p n + B X B 2 > 0, for n > , (41) 
AxA 2 p 2n + 2p 2n+1 + 4A 2 p n+1 +2p+ [(1 + 4n)Ax + (1 - 4n)A 2 A] p n + A > 0, for n < . (42) 

We address here the evident conditions for which the WEC is satisfied as follows: 

* Bx > 0, B 2 > for 0<n<l/4, 

* Ax > 0, A 2 > 0, for -l/4<n<0. 
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It is obvious that these conditions are not unique. For n > (n < 0), the necessity of plotting 
the function Ap 2n + 2p 2n+1 + AB 2 p n+1 + 2B|p + [(1 + An)AB 2 + (1 - 4n)S x ] p n + B X B 2 , (A 1 A 2 p 2n + 
2A\p 2n+1 + 4A 2 p n+1 + 2p + [(1 + 4n)Ai + (1 - An)A 2 k] p n + A) varying two of the input parameters 
and fixing the other, appears when at least one of the following expressions is negative, (B 2 , B\B 2 or 
(1 + 4n)AB 2 + (1 - 4n)Bi), (AiA 2 , A 2 , or (1 + 4n)Ai + (1 - 4n)A 2 A). We present some examples of 
theses cases in Fig. 2. 

• The SEC 

The strong energy condition is realized by combining the NEC with p e // + %Peff > 0. This latter 
reads, 

Peff+3 Peff =p + 2pf p -2f(p)>0 . (43) 

Making use of the expressions in (J6)), one obtains a fraction whose the denominator is always positive 
and the numerator reads 

p 2n+i + 2B 2 p n+1 +B 2 p-2 [(1 - n)AB 2 + (1 + n)B x ] p n - 2Ap 2n - 2B X B 2 > , for n > , (44) 
A 2 p 2n+1 + 2A 2 p n+1 + p - 2 [(1 - n)Ai + (1 + n)A 2 A] p n - 2A 1 A 2 p 2n - 2A > for n < . (45) 

Now, combining (|4~4")) - (l4"5"j) with the NEC, on gets the following conditions for the SEC 

2p 2n+1 - 2Ap 2n + AB 2 p n+1 + 2B\p - 2 [(1 - 2n)AB 2 + (1 + 2n)B{\ p n - 2B Y B 2 > , for n > , (46) 
2Alp 2 " +1 - 2AiA 2j o 2 ™ + 4A 2/ 3 ,l+1 + 2p - 2 [(1 - 2n)A 1 + (1 + 2n)A 2 A] - 2A > , for n < . (47) 

In this case, there is any obvious condition for satisfying the SEC. However, values can be found, by 
plotting 2p 2n+1 - 2Ap 2n + AB 2 p n+1 + 2B\p - 2 [(1 - 2n)AB 2 + (1 + 2n)B Y \ p n - 2B Y B 2 and 2A 2 2 p 2n+1 - 
2AxA 2 p 2n + AA 2 p n+1 + 2p - 2 [(1 - 2n)A 1 + (1 + 2n)A 2 A] p n - 2A in terms of two of the parameters. 
Some examples for illustrating some of these cases are presented in Fig. 3. 

• The DEC 

The dominant energy condition is characterized by the WEC combined with p e ff—p e f / > 0. Following 
the same steps as in the previous cases, one easily obtains the DEC as 

p 2n+i + Ap 2n + 2 B 2 p n+1 + B\ p + B X B 2 + \{2n + 1) AB 2 + (1 - 2n)Bi] p n > for n > , (48) 
A 2 p 2n+i + A^p 211 + 2A 2 p n+1 +p + A+ [(2n + l)A 1 + (1 - 2n)AA 2 ] p n > for n < . (49) 

The evident conditions read 

*B 1 >0, B 2 >0 and n > 1/2 . 

Evidently, other conditions may lead to the accomplishment of the DEC, but, only plotting the 
functions in (|4"8]) - (|49"|) . This can occur when at least one of the following expressions is negative, B 2 , 
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£?i£?2, (2n + l)Ai?2 + (1 — 2n)B\ for (|48|> . We recall that the negative parameters must be suitably 
chosen without loosing the positivity of the expressions (|48 |) -(|49 |) . We present some of theses cases in 
Fig. 4. 



4.2 Studying the case f(T) = a 3 In 9 (b 3 T m ) 

Here we will work with the fundamental conditions for which the model allows the avoidance of the 
Big Rip. So, we propose to check if the range of parameters for which the singularity may be cured can 
also make the model satisfying the energy conditions. Here, the first derivative of f(p) also plays an 
important role. Deriving /(p) with respect to the energy density p, one gets 

fM = ^ h £- l (b 3 p™) . (50) 

We believe that each step of constructing the four energy conditions is now clear and we simply present 
the results and comments as follows: 

• NEC 

p + 2qma 3 + p\n q - 1 (b 3 p m )>0 . (51) 

The evident conditions for obtaining this are 

• q > 0, to < 0, a 3 < 0, with p m > l/b 3 . 

It is important to note that this list is not exhaustive, since in other conditions different from the above 
ones, the NEC could still be realized. This is a set of the situations where for example 2qma 3 < 0, but its 
absolute value is less than the absolute value of p+p\n q ~ 1 (b 3 p m ), i.e., |2gma3| < Ip+pln 9 " 1 (b 3 p m )\. This 
situation requires knowing some intervals to which the parameters must belong. We present this feature 
by plotting the function corresponding to the expression (I5ip in terms of some of the input parameters 
fixing the other. See the graph at left side in Fig. 5. 

• WEC 

2p + 2qma 3 + (p + 2qma 3 ) In 9 " 1 (b 3 p m ) + a 3 In 9 (b 3 p m ) > . (52) 

In this case by plotting the function ([52"T) . the WEC can be realized graphically. This is the set of 
situations where one of the term in the sum (f52"j) is negative, but it absolute value is less that the 
absolute value of the sum of the other. An example for this is when 2qma 3 < 0, but |2gma3| < 2p + 
| (p + 2qma 3 ) In 9-1 (b 3 p m ) + a 3 In 9 (b 3 p m )\. See the graph at right side in Fig. 5. 

• SEC 

2p + 2qma 3 + {p + 2qma 3 ) In 9 " 1 {b 3 p m ) - 2a 3 In 9 (b 3 p m ) > . (53) 

In this case, evident constraints on the input parameters in order to realize this energy conditions are 
presented as follows: 
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* q > 0, m < 0, a 3 < 0, with p m > l/b 3 , 

As presented in the previous cases, other conditions may also realize this energy conditions. This can 
be observed by plotting the function in (|53[) in terms of some input parameters, fixing the other. See the 
graph at left side in Fig. 6. 

• DEC 

■ip + 2qma 3 + (p + iqma 3 ) In 9-1 (b 3 p m ) + 3a 3 In 9 (b 3 p rn ) > . (54) 

Here, constraints may also lead to the DEC, but this is clear by plotting the function ([54]) . as in the 
previous cases. We present an illustrative example at right side in Fig. 6. 

We mention that for all the graphs, the parameters are normalised to 10 -121 Planck units. Remark 
that the current value of the cosmological constant is about 1.7 x 10 -121 and the energy density of the 
usual matter is about p — 0.1 x 10~ 121 [Tl]. Then, with the normalization, we get A = 1.7 and p = 0.1 
for the cosmological constant and the energy density of the usual matter respectively, which are the 
values used for plotting the graph in the figures. Moreover, we observe that all the conditions which lead 
to the accomplishment of the NEC also satisfy the WEC, the SEC and the DEC. This is an obvious 
deduction from (|29|) - ([32|) . In the same way, it is easy to observe that all the conditions for which the 
WEC is satisfied, also lead to the accomplishment of the DEC. From the same expressions (f3"0"]) -(f3"T j) . we 
see that the sufficient condition for the requirement of satisfying the WEC leads to the accomplishment 
of the SEC is f(T) < 0. Thus, among the conditions for which the WEC is satisfied, those for which the 
function is always negative, also lead to the accomplishment of the SEC. However, about the SEC with 
respect to the DEC, the situation is quite different. The requirement for, the conditions which satisfy the 
SEC, leading to the accomplishment of the DEC, is always having f(T) > 0. This means that, among the 
conditions for which the SEC is satisfied, those for which the function f(T) is always positive, guarantee 
the accomplishment of the DEC. 



5 Perturbations and stabilities in R + 2/(T) gravity 

In this section we propose to study the perturbations around the models used in this work. We can 
start establishing the perturbed equations for the case R+ 2/(T), but the two models will be studied as 
specific cases. 

For this purpose, let us assume a general solution for the cosmological background of FRW metric, 
which is given by a Hubble parameter H — H^t) that satisfies the background equation (|27|) using (|25|) . 
for R + 2/(T) gravity. The evolution of the matter energy density can be expressed in terms of this 
particular solution by solving the continuity equation around Hb(t), 

p b + 3H b (t)p b = , (55) 



13 



yielding 

p b {t) = p Q e-^ H ^ . (56) 

We recall that we are considering that the ordinary content of the universe is pressureless. Since we 
are interesting in studying the perturbations around the solutions H = H b (t), we will consider small 
deviations from the Hubble parameter and the energy density, i.e., we can write the Hubble parameter 
and the ordinary energy density as |23j 

H(t) = H b (t)(l + 6(t)) , P (t) = Pb (t)(l + S m (t)) . (57) 

In order to study the behaviour of these perturbations in the linear regime, we expand the function f(T) 
in powers of T b (or p b ) evaluated at the solution H = H b (t), as 

f(T) = f( P ) = f b + f b p (p-p b )+0 2 , (58) 

where the superscript b refers to the background values of f(T) and its derivatives evaluated at T = T b 
(or p = p b ). Here, the O term includes all the terms proportional to the square or higher powers of T (or 
p). Then, only the linear terms of the induced perturbations will be considered. Hence, by making use 
of the expression (|57|) in the Eqs. (|27|) and (l25j) . one gets the equation for the perturbation 5{t) in the 
linear approximation, 

(Pb + 3p b f b p + 2p\f pp ) 5 m (t) = 6HlS(t) . (59) 
On the other hand, there is a second perturbed equation from the matter continuity equation, 

5 m + SH b {t)S(t) = . (60) 
By combining Eqs. (|59")l and (f6"U| one gets the following equation for the matter perturbation 

2H b S rn + (p b + 3p b f b p + 2p 2 b f b pp ) S m = , (61) 

from which we obtain 

S m (t) = Ci exp I - i J C b dt} , C b = -g- (1 + 3f p + 2 Pb f b pp ) , (62) 
where Ci is an integration constant. By using the relation (|60p . the perturbation 5 reads 

^*) = -§^ ex p|-^ / C bdt \ ■ (63) 



6H b r [ 2, 

Let us now consider two cosmological solutions and analyse their stability by the use of the models treated 
in this work: de Sitter solutions and power-law solutions. 
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5.1 Stability of de Sitter solutions 

In de Sitter solutions, the Hubble parameter is constant and one has 

H b (t)=H ,^a(t)=a Q e Hot , 



(64) 



where Hq is constant. 

With this scale factor, the energy density of the background becomes p b — poe~ 3Hot , with which one 
has dpi, — —3Hopi,dt. By using this, one can cast the integral in (|62p into 

1 



C h dt = -- 



3H 2 



(1 + 3/ p 6 + 2p b f b pp ) dp b 



(65) 



AT"+Bi 



• Treating the model f(T) — T „ +B2 . 

This case corresponds to n > 0, and the integral (|6"5"|) can be expressed as 



C b dt 



and Ch is written as 



3H 2 l 



Pb 



AB 2 -B 1 2n(AB 2 -B 1 )p% 



Pb 



B 2 + P n b (B 2 + p^Y 

n{2n + 1)(AB 2 - B 1 )p% 4n 2 (AB 2 - B x )p, 



n _ n -3nH t 

P b — Po e 



{P b +B2? 



(66) 



(67) 



(P b + B 2 f 

We see from (|66|) and (l67|) that for n > 0, and as the time evolves, the stability of de Sitter solutions 

aiT n + bl 
a 2 T"+fc 2 : 



requires B 2 ^ 0. In other word, for the initial model f(T) — ^^sxfe , de Sitter solutions are stables if 



and only if a 2 =^ and b 2 =/= 0. 

• Treating the model /(T) = A ^ n %\ ■ 

This case corresponds to n < 0, and the integral (|6"5|) . multiplied by —1/2, can be expressed as 



I I C b dt = 



1 



6H 2 



Pb 



A\ - AA 2 2n(A 1 - AA 2 )p[ l 



, P b = Po e 



^ — 37lHnt 



and Cb is written as 



Cb = w 



Pb 



A 2 {l + A 2P n b ) (l + A 2 p%r 
n(2n + \){A X - AA 2 )p% ±n 2 A 2 (A X - AA 2 )p 2 b n 



(69) 



{A 2 pl + l) 2 {A 2 pl + 1)3 

Here, for n < 0, as the time evolves, both (|6"5|) and (j6"9")l tend to +oo. Thus the perturbation will grow 
exponentially, and this particular de Sitter solution becomes unstable. Note that this result does not 
depend on any of the parameters A\ or A 2 . 

• Treating the model f(T) = a 3 \n q (b 3 T m ). 

With this model, the integral (|6"5j) . multiplied by —1/2, can be performed and one gets 



\ I ° bdt = 6B~o^ b + ° 3 + 2mqa3 lnQ l {b3P " 1 ^ ' 



with the corresponding expression of C b being 
1 



[p b + qma 3 In 9 " 1 (b 3 p m ) + 2mq(q - l)a 3 In"" 2 (b 3 p m )] 



(70) 



(71) 



15 



Let us recall that this model (f(T) — 03 In 9 (bsT" 1 )), leads to the avoidance of the Big Rip for q > 
2a/ {a — 1) and m < 0, where a > 1, as we have previously shown. These conditions also allow the 
model to satisfy the energy conditions. Now, let us check what happens about the stability with these 
conditions. First, note that the relation q > 2a/ (a — 1) can be cast into q > 2 + 2/(a — 1), showing that 
q > 2 because of a > 1. By choosing 0,3 < 0, we see that, within the conditions q > 2 and m < 0, the 
expressions (170)) and (I7T|) tend to —00 as the time evolves, and this ensures the decay of the perturbation, 
leading to the stability of de Sitter solutions with this model. Thus, regarding to the stability of dc Sitter 
solutions, the energy conditions and the late time acceleration, provided with the conditions 03 < 0, 
&3 > 0, q > 2 and m < 0, we can conclude that the model may be cosmologically acceptable. 



5.2 Stability of power-law solutions 



As we are dealing with dust as ordinary content of the universe, we will be interested to the scale 
factor 



a{t)a t 2/3 -> H b (t) = ^ , Pb oc t~ 2 

• Treating the model f(T) = ^C+g 1 . 

In this case, n > 0, and one can perform the integral 



- \ fc b dt = -^- 

2 J 4a 



. ... (3n + 4)(AB 2 - B 1 )t 2 (AB 2 - Bi)t 2 
In (t) h 



B 2 (l + B 2 t 2n ) 



(4n + 7)(Ag 2 -g 1 )t 2 (I 
+ 2B 2 2J Hn' M 



2B 2 n{\ + B 2 t 2n 

-;-B 2 t 2n 
n 



with, 



6 2a 



1 2n(n+l){kB 2 - Bx)r 



4n 2 (AB 2 - Bx)r 



t {t-^ + B 2 f {t-^ + B 2 f 

where we have set p$ = 1, and 2 F\ is the hypergeometric function defined by 

(Al) r (A 2 ) r Z r 



2 Fi(Ai,A 2 ,A 3 ;z) = ^ 



r=0 



(A 3 )r 



(A) r = A(A + 1)(A + 2). ..(A + r - 1) , (A) = 1 



(72) 



(73) 
(74) 

(75) 
(76) 



As the time evolves, conditions are required for guaranteeing the decay of the perturbation. For B 2 < 0, it 
is necessary to have Ai3 2 < Bi, which means that B\ can be positive, or negative but with |AS 2 | > \B\\. 
In the case where B 2 > one may observe two sub-cases, i.e., for an even r, and an odd r . For an 
even r, as the time evolves, the necessary condition for guaranteeing the decay of the perturbation is 
A.B 2 > 61, meaning that the parameter B\ can be negative, or positive. On the other hand, for an odd 
r, the requirement for getting the decay of the perturbation is A_B 2 < £?i, meaning that B\ > 0. 
• Treating the model f(T) = ^ffjpf ■ 
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Here, n < 0, and the integral can be performed as 

(2n + 3)(A 1 - AA 2 )t 2 n(A 1 - AA 2 )t 2 



r / C h dt = 



3 
4a 



ln(<)- 



2{A 2 + < 2 «) 



+ 



?,t 2 {A 1 - AA 2 
2A 2 



■ 2 -Pi 



(A 2 +i 2 «) 2 
1 t 2 ™ 



with 



C h = 



2a 



n(2n + - AA 2 )^ 2n+1 4A 2 n 2 (A 1 - AA 2 )t 



-4n+l 



(77) 



(78) 



(A 2 t- 2n + l) 2 (A 2 t- 2 » + l) 3 

As the time evolves, the argument of the hypergeometric function tends to zero and the hypergeometric 
function tends to 1. Thus, the dominant term in (|77[) reads 

3 



:(A 2 -n)(Ai-AA 2 y 



(79) 



Here, one can distinguish two cases: (A 2 — n > and A\ < AA 2 ) and (A 2 — n < and A\ > AA 2 ). In 
the first case, one gets A 2 > n meaning that A 2 can be positive, or negative but with \A 2 \ < \n\. When 
A 2 > 0, Ai can be positive or negative, due to the relation Ai < AA 2 , while for A 2 < 0, A\ is necessarily 
negative. In the second case, one gets A 2 < n, meaning that A 2 < 0, which allows Ai to be positive, due 
to the relation A\> AA 2 . 

We observe that some of the conditions for which the stability occurs, are also compatible with some 
energy conditions. This shows that for some values of the input parameters, acceptable models can 
be obtained, at least regarding to the energy conditions, the stability, the late time acceleration of the 
universe and the avoidance of the Big Rip. 

• Treating the model f(T) = a 3 In 9 (b 3 T m ). 

As we have done in the previous cases, the integral can be performed, yielding 



-\ic h = - 

with 



4a 



In (t) + l-a 3 qbf m^[q, i®] + a 3 q(q - l)m«bj T[q - 1, ^] 
2 m TYi 



3 
2a 



+ a 3 qm t In 9 " 1 (b 3 t~ 2m ) + 2a 3 q(q - l)m 2 t In 9 " 2 (6 3 r 



■2m \ 



? (i)=ln(M- 2m ) -(80) 



(81) 



As we have previously shown, this model cures the Big Rip for q > 2 and m < 0. With these conditions, 
as the time evolves, only the term In (t) grows. Since —3 In (t) / (4<zo) is negative for large value of the time, 
it is easy to observe that the perturbation decays, and this corresponds to the stability of the power-law 
solutions with this model. Observe that in this case, the constraints on the parameters q and m for which 
all the energy conditions are satisfied, leads to the stability of the power-law solutions. Thus, regarding 
to the stability, the energy conditions, the late time acceleration of the universe and the avoidance of the 
Big Rip, we can conclude that this model can be cosmologically acceptable for a 3 ^ 0, b 3 > 0, q > 2 and 
m < 0. 
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6 Discussions 



In this work we discussed the viability of an interesting alternative gravitational theory, namely, 
modified f(R,T) theory of gravity, where R is the curvature scalar and T the trace of the energy mo- 
mentum tensor. For our purpose we took the case where the algebraic function f(R, T) is cast into 
f(R,T) = fx(R) + / 2 (T), and focused our attention to the assumption f(R) = R and f 2 (T) = 2/(T). 
The aim of the work is to analyse the viability of two models, f(T) — (aiT" + 6i)/(aaT™ + b 2 ) and 
/(T) = 03 In 9 (b^T" 1 ), (where 01, &2, b±, 62, ti, 0,3, 63, q and m are input parameters) by investigating 
for which range of each parameters, the models satisfy all the energy conditions or some of them. We 
also assumed that the usual matter content of the universe is pressureless (dust), for which the trace 
T = p. Depending of having 02 7^ or 62 ^ 0, and according to the fact that the present stage of the 
universe may be characterized by the ACDM model, we distinguished two sub-models for the first one. 
Thus for a2 ^ and n > 0, the first model reduces into f(p) — (Ap n + B\)/{p n + B2), where B\ = &i/a 2 
and B 2 = ^2/^2, while for b 2 7^ and n < 0, the model reduces into f(p) = (Aip n + A)/(^4 2 p" + A), 
where A\ — a\/b% and A 2 — 0,2/^2 ■ Energy conditions are generalised to the f(R,T) by acting directly 
on the effective energy density and pressure, respectively p e ff and p e ff of the corresponding dark fluid. 
The same energy conditions, the null, weak, strong and dominant, applied to the ordinary energy density 
and pressure in GR are then valid for the effective energy density and pressure. Besides to the evident 
constraints on the input parameters for each of the expressions of f(T), for obtaining R + 2/(T) gravity 
models that satisfy energy conditions, we also presented figures characterizing models satisfying energy 
conditions, and whose parameters cannot be determined analytically, but rather, numerically. For plot- 
ting the graph, the current value of the cosmological constant is used A ~ 1.7 x 10~ 121 Planck units [T3], 
and the value of the energy density of the usual matter is considered as p ~ 0.1 x 10 -121 Planck units. 
The input parameters are normalized to 10~ 121 Plank units, such the cosmological constant and the 
energy density used for plotting the graph are 1.7 and 0.1, respectively. We also discussed the situations 
where all or some of the energy conditions are satisfied. We observe that all the conditions which lead 
to the accomplishment of the NEC also satisfy the WEC, the SEC and the DEC. This is an obvious 
deduction from (|29 |) - ([32"T) . In the same way, it is easy to observe that all the conditions for which the 
WEC is satisfied, also lead to the accomplishment of the DEC. From the same expressions ff50|) - (f5Tj) . we 
see that the sufficient condition for the requirement of satisfying the WEC leads to the accomplishment 
of the SEC is f(T) < 0. Thus, among the conditions for which the WEC is satisfied, those for which the 
function is always negative, also lead to the accomplishment of the SEC. However, about the SEC with 
respect to the DEC, the situation is quite different. The requirement for the conditions that satisfy the 
SEC leading to the accomplishment of the DEC, is always having f(T) > 0. This means that, among the 
conditions for which the SEC is satisfied, those for which the function f(T) is always positive, guarantee 
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the accomplishment of the DEC 

Moreover, in order to check the viability of the models in some ways, we investigate the stability of 
de Sitter solutions and power-law solutions under perturbations in the FRW metric in the framework 
of R + 2/(T) gravity, by using the models proposed in this work. The results show that only for the 
models f(T) = (AT" + Bi)/(T" + B 2 ) and f(T) = a 3 In 9 (b 3 T m ) de Sitter solutions can present stability. 
However, for the power-law solutions, the stability can be obtained for each model under some conditions. 
We observe that some of the conditions that lead to the stabilities of de Sitter solutions and power-law 
solutions are compatible with the accomplishment of some of the energy conditions, the late-time cosmic 
acceleration and the avoidance the type-I future finite time singularity (Big Rip). We conclude that, in 
the framework of R + 2f(T), viable models of type f(T) = l^tbl and /( T ) = a 3 ln<? ( 6 3 T ™) can be 
obtained, when the usual matter is assumed as pressureless. 

Despite to these satisfactory results obtained in the framework of R + 2f(T) gravity, it is important to 
recall that this work may be view as a particular case of the general one where f(R, T) = fi(R) + f2(T), 
with fi(R) 7^ R. Our purpose for a coming work is to undertake this general case, for which the f(R) 
gravity is also a particular case, and analyse the energy conditions, the stability of some known scale 
factor solutions, against homogeneous perturbations and also the cosmological background evolution. 
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Figure 1: The graph at the left side of this figure corresponds to the NEC by plotting 2np n_1 (AB 2 - B\) + (p n + B 2 ) 2 
in functions of B\ and B 2 with p = 0.1, n = 1, A = 1.7. This case is different from the evident conditions. Here, 
2np n ~ 1 (AB2 — B\), but the NEC is realized. The graph at the right side of the figure corresponds to the NEC but in this 
case by plotting 2np n ~ 1 {A\ - A 2 A) + (A 2 p n + l) 2 in functions of A\ and A2 with n = -1, A = 1.7 an p = 0.1. This 
situation is also different from those of the evident conditions, but also realizes the NEC. 




Figure 2: The graph at the left side is obtained by plotting the function Ap 2 ™ + 2p 2n+1 + AB 2 p n+1 + 2B 2 p + 
[(1 + 4n)AB 2 + (1 - 4n)Bi] p n + B\B 2 in terms of suitable values of B\ and B 2 , with n = 1, A = 1.7 and p = 0.1. 
We have a situation different from the evident ones but the WEC is satisfied. At the right side of the figure one has the 
graph obtained by plotting Ai A 2 p 2n + 2p 2n+1 + 4A 2 p n+1 + 2p + [(1 + 4n)Ai + (1 - 4n)A 2 A] p n + A in terms of Ai and 
A 2 , with n = —1, A = 1.7 and p = 0.1. We also have here a case different from the evident ones but for which the WEC is 
still realized. 
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Figure 3: At the left side the SEC is represented by plotting 2p 2n+1 - 2Ap 2n + AB 2 p n+1 + 2B\p - 
2 [(1 - 2n)AB 2 + (1 + 2n)B 1 ] p n - 2B\B 2 in functions of Bi and B 2 , setting A = 1.7, p = 0.1 and n = 1. The graph at right 
side of the figure represents the SEC where 2A 2 p 2n+1 - 2Ai A 2 p 2n + &A 2 p n+1 + 2p - 2 [(1 - 2n)A 1 + (1 + 2n)A 2 A] p n - 2A 
is plotted in terms of Ai and A 2 , setting n = —1, p = 0.1 and A = 1.7. 




Figure 4: The graph at the left side of this figure characterizes the DEC, coming from p 2n + 1 + Ap 2 ' 1 + 2B 2 p' I+1 + 
B\p + BxB 2 + [(2n + 1)A_B 2 + (1 - 2n)B 1 ] p" plotted in functions of B 1 and B 2 , by setting n = 1, p = 0.1 and A = 
1.7, while the graph at the right side represent the DEC from the plot of A?,p 2n+1 + A\A 2 p 2ri + 2A 2 p n+1 + p + A + 
[(2n + l)Ai + (1 - 2n)AA 2 ] p" in functions of Ai and A 2 , with n = -1, p = 0.1 and A = 1.7. 
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Figure 5: The graph at the left side of this figure represents the NEC and is obtained by plotting p+2qmas +p In 9-1 (&3p m ) 
in functions of m and 63, setting p = 0.1, 03 = 1 and q = 3. This case is different from the evident conditions but still realizes 
the NEC. The graph at the right side of the figure corresponds to the WEC, plotting 2p+2qma-j + (p + 2qma^) In 9-1 (&3P m ) + 
CI3 In 5 (fe3P m ) in functions of m and 63, where we set q = 3, 03 = 1 and p = 0.1. 




Figure 6: The graph at the left side characterizes the SEC is obtained by plotting 2p+2qma^ + (p + 2qma;i)ln q ~ 1 (63P" 1 ) — 
203 ln q (b3p m ) in functions of m and £>3 where we set p = 0.1, as = 10 — 3 and q = 3. This situation is different from the 
evident conditions, but also realizes the SEC. The graph at the right side of the figure is obtained by plotting 3p + 2qmas + 
(p + iqma-j) In 9-1 (bsp m ) + 3123 In 9 (bsp m ) in functions of m and 63, with q = 3, as = 1 and p = 0.1. This graph represents 
the DEC. 
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